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Analytic computations of massive one-loop amplitudes 
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We show some new applications of on-shell methods to calculate compact helicity amplitudes for tt production 
through gluon fusion. The rational and mass renormalisation contributions are extracted from two independent 
Feynman diagram based approaches. 



1. Introduction 

With the first data from the LHC now arriving 
the need for full NLO studies of multi-leg QCD 
processes is becoming ever more important. The- 
oretical developments in recent years have rapidly 
opened the range of processes that can be feasibly 
considered. Modern methods using combinations 
of unitarity, factorisation and algebraic tensor re- 
duction have been able to produce automated 
numerical codes for phenomenological studies of 
2 — >■ 4 processes jlj-iS]- A more detailed overview 
of the current status of NLO studies can be found 
inref. ^. 

One of the most important areas for both phe- 
nomenological applications and as a benchmark 
for theoretical complexity is the inclusion of full 
mass dependence in NLO QCD studies. Numeri- 
cal methods have been at the forefront of progress 
in this direction and impressive results have al- 
ready been achieved for pp ^ tt + j 0, [Hj , 
pp ^ tt + bb d, H, [13] and more recently pp — 
tt + 2j f^. 

Progress in fully analytic computations of mas- 
sive one-loop amplitudes has been relatively slow. 



Despite many theoretical developments |13l - ll7| . 
there are relatively few complete calculations 
used in phenomenological studies. The hope that 
analytical computations would provide faster and 
numerically more stable results motivates us to 
study this problem in more detail and investigate 
to what extent the problem can be automated. 
In this contribution we present new compact an- 
alytic one-loop amplitudes for gg — )■ tt scatter- 
ing. In section [2] we briefly outline the on-shell 
techniques we have used and discuss the problem 
of mass-renormalisation within this approach. In 
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section [3] we discuss two independent Feynman 
calculations which have been used to verify our 
results and generate compact representations of 
the rational terms. We define our conventions for 
the Spinor- Helicity formalism in section |4] before 
presenting our results in[5j We discuss some pos- 
sibilities for alternative methods for determining 
tadpole coefficients before reaching our conclu- 
sions. 

2. Generalised Unitarity 

Generalised unitarity has become an essential 
tool in the computation of multi-particle one-loop 
amplitudes over the last few years. Building from 
the pioneering work of Bern, Dixon, Dunbar and 
Kosowcr in the mid-nineties fisf the modern pic- 
ture of unitarity gives a purely algebraic approach 
to the computations of one-loop amplitudes mak- 
ing use of complex analysis |19H23| . 

Quadruple cuts with complex momenta com- 
pletely freeze the four dimensional loop integra- 
tion leading directly to the algebraic evaluation of 
the associated scalar box integral coefficient (l9j . 
The technique of integrand reduction first pro- 
posed by Ossola, Papadopoulos and Pittau [l^l 
was shown to apply elegantly to the extraction 
of triangle and bubble coefficients using Laurent 
expansions of the unconstrained integrations by 
Forde [2lj. This method has been generalised to 
the case of massive amplitudes [13] and is the ba- 
sis for the method used here. 

In these proceedings we apply the technique 
of generalised unitarity to compute compact an- 
alytic expression for one loop tigg helicity am- 
plitudes. The results are obtained from a semi- 
automated system using FORM 24] and Maple. 
The integral basis and integrands are generated 
with Maple before the tree amplitudes (from 
BCFW recursion relations 2^ 2^) are expanded 
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using FORM to extract closed analytic forms for 
the integral coefficients. IR consistency equations 
have been verified and used to constrain the form 
of the amplitude and ensure a compact represen- 
tation. 

The colour ordered amplitudes can be ex- 
pressed in Z? = 4 — 2e dimensions as a product of 
rational coefficients multiplied by divergent scalar 
integrals, 



K3 



(1) 



In the massive amplitudes the wave-function 
renormalisation and tadpole contributions are 
non-zero and require mass-renormalisation to 
be performed before the terms proportional to 
log(m^) can be determined. This presents a prob- 
lem when formulating a purely on-shell construc- 
tion of the amplitude since cuts in the log(m^) 
channels contain explicit divergences. One can 
circumvent this problem by explicitly removing 
the singularities but this causes gauge invariance 
to be broken in intermediate stages of the com- 
putation 27[. In 4 — 2e dimensions the univer- 
sal pole structure is sufficient to determine the 
log(r7i^) terms however additional rational terms 
are introduced, 



K' 



n + o{e). 



(2) 



where the sum over K2 includes only the "off- 
shell" bubbles with well defined unitarity cuts. In 
our approach both C2;m^ and TZ are computed us- 
ing Feynman diagram based techniques. We can 
then verify that the coefficient of C2-m^ matches 
that predicted by the IR constraints. We then 
demonstrate for the example of the t~^g~^g^t~^ 
amplitude that R^^ can also be determined us- 
ing the cut-constructiblc information. 



3. Feynman diagram approach 

The amplitude was computed using two inde- 
pendent approaches. An analytic computation 
using traditional Passarino-Veltman tensor re- 



duction was performed with the help of DIANA [28 1 



(diagram generation) and FORM [2j|. A second 
method employed a numerical implementation of 
tensor reduction based on recurrence relations in 

shifted dimensions [ii-Iil. 



The 37 colour-ordered one-loop diagrams con- 
tributing to the process gg — >■ tt, including ghosts 
and on-shell mass renormalisation counter terms, 
were further processed analytically in two inde- 
pendent FORM codes to generate tensor integral 
representations of the amplitudes. The calcula- 
tions were performed using the spinor-helicity for- 
malism in the Four-Dimensional Ilelicity(FDII) 
regularisation scheme. 

Full agreement between the unitarity and Feyn- 
man approaches was achieved with the use of two 
independent CH — I- implementations of the spinor 
products and the qcdloop package for finite parts 
of the 4 — 2e-dimensional scalar integrals 133| . 

4. Spinor/Helicity Formalism 

For massless particles it is possible to com- 
pletely decompose all momenta into a basis of 
two-component Weyl spinors since, 



= \{p\r\p]- 



(3) 



The polarisation vectors and fermion wave- 
functions then fit easily into a helicity basis. 



u+ip) = \p) 



V2(Cp) 



u^ip) - \p] (4) 

{p\rm 



(5) 



Kleiss and Stirling described how to construct 
well defined helicity states for massive momenta 
through the introduction of an arbitrary massless 
vector which defines the reference frame Q| . For 
a massive vector P we use a massless vector rjp 
to define a massless projected vector P'': 



2P-r] 



(6) 
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The u and v spinors can then be defined by: 



u±{P,m;P\'qp) = — . „, , (7) 



«±(F,m;F^^7p) 



(77P T\P'±) 
{f-m)\r]p±) 



(8) 



The freedom to keep rjp arbitrary gives us the 
ability to relate the positive and negative helicity 
states by applying a simple transformation [P^ ^ 
tip)- 



v_{P,m;P\r,p) = 



v+{P,m-r^p,P^). (9) 



It is therefore sufficient to consider helicity am- 
plitudes where all massive fermion have positive 
helicity and obtain the remaining amplitudes by 
applying eq.®. For the four-point ttgg exam- 
ple considered here this gives us two independent 
configurations, + + +-f- and + H h. 



4;i(l„2,3,4f) = ^m(l„2,3,4f) 
- -L^M(1,,2,3,%- ^A[/l(l„2,3,4f) 

-^A[«l(l*,2,3,4,), 

4^](lt,4f;2,3)= {^'^'(1^,2,3,4,-) 

P(2,3) ^ 



AW(lt,2,4f,3) + A[«l(lt,2,3,4f 



(12) 



For the sake of brevity we restrict ourselves to the 
+ + ++ helicity configuration. The remaining in- 
dependent helicity amplitude has been generated 
and checked against the literature but will be pre- 
sented elsewhere. The tree level amplitude can be 
written as, 



4'')(l+,2+,3+,4+): 
. 3 [23] (771 774) 



(23)(2|l|2](ryil^)(r;44^)- 



(13) 



5. ttgg helicity amplitudes 

In this section we demonstrate the techniques 
by applying them to the well studied case of top 
pair production through gluon fusion. The final 
results are obtained in an extremely compact an- 
alytic form and are in agreement with the results 
in the literature sB" 33| • The amplitude is decom- 
posed into colour ordered primitive amplitudes as 
outlined in 13811: 



We find it convenient to write down the one-loop 
amplitudes in terms of the scalar integrals (us- 
ing the notation of 33]) together with two finite 
integral functions and F2- By defining, 

F4 (0, TO^, m^, 0, S12, S23i m'^,m^, 0, m^) = 
li (0, TO^, m^, 0, S12, S23, m'^,m'^, 0, to^) 



(2|1|2] 



(14) 



4"^(lt,2,3,4f) = 
J2 (r'^^T'^^)nu4°\lt,2,3,4,-) (10) 

P(2,3) 

4^^(lt,2,3,4f) = 

E A^(r"=r"-^)»:u4a(it,2,3,4f) 

P(2,3) 

+ ^'^^'^M,,,,4;](l„4,-;2,3), (11) 



we absorb all the divergent ilog(^^), with 

/3 — ~ ^I^' ^-'^^ three-mass triangle coef- 
ficient. We can then make use of the IR constraint 
on the integral coefficients to simplify the result: 

*^4a|2|3|4 + (2|1|2]C^^|3|4;^ = 

(s23-2m2)(2|l|2]A(")(li,2,3,4f). (15) 
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We also move the log(m^) dependence into the 
on-shell bubble contributions by defining, 

F2 (si2,0,m2) = 

h (s23,m^,m^) - I2 (0, 



[27| using: 



m^m2). (17) 



The remaining inte gral s have been computed an- 
alytically long ago (39I |40| . We have made use 
of qcdloop '335 in making numerical comparisons 
with the literature. The analytic expressions for 
the leading colour primitive amplitudes are given 
in equations eqs. (|18ll2n|) . 

6. Tadpole coefRcients and IR/UV con- 
straints 

The universal IR and UV behaviour provide 
strong constraints on the form of the terms pro- 
portional to log(TO^). It has been proposed long 
ago by Bern and Morgan that this information 
can be used to completely fix the on-shell bubble 
and tadpole coefficients 13j- Here we show, for 
the + + H — h example considered above, how this 
can be extended to deal with massive external 
particles. It is hoped that proceeding along simi- 
lar lines would give a general method for any mas- 
sive amplitude. Using the universal pole structure 



4l|,i42] the overall coefficient of 1/e — log(m ) can 



be described by: 

AW(lt,2,3,4f) = 
A«'-(lt,2,3,4,-)-f 7^(lt,2,3,4,-) 
Cf ( Ai| 



2,3,4,) — 



(1 - 2e) \m 



(21) 



In our example we find: 



C2|„2(l+,2+,3+,4t) = 
-2CfC2;12(1^2+,3+,4+) 
Ci(l+ 2+,3+,4±) = 
-Cf^AW(l+ 2+,3+,4+). 



(22) 
(23) 



R 



(l+,2+,3+,4t 



7e(l+, 2+, 3+, 4+) + -A(0)(1+, 2+, 3+, 4+), 

(24) 



where 7?. is defined according to eq. The 
equation applies to both Left- and Right-moving 
primitive amplitudes. We note that this simple 
relation holds only for the -|- -I- helicity ampli- 
tude since the "on-shell" bubble happens to can- 
cel the "off-shell" bubble precisely. For the gen- 
eral case a more complicated relation would need 
to be determined. 

7. Outlook 

We have shown that combinations of gener- 
alised unitarity and optimised Feynman diagram 
computations can be an efficient method to gen- 
erate compact analytic expressions for massive 
one-loop amplitudes. We have shown complete 
helicity amplitudes for tt production via gluon 
fusion in the "all-plus" configuration. We have 
checked our results against those in the literature 
and found complete agreement in all cases. We 
hope that the result presented here will allow for 
a faster evaluation of the NLO cross-section. 

The methods employed are quite general and 
have also been successful in studying the more 
complicated helicity configurations. However the 
full results for this process, including the sub- 
leading colour contributions, will be presented in 
a forthcoming publicatior0. 
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One can then find agreement with the D- 
dimensional rational term as defined in reference 



^The full amplitudes have been subsequently presented in 
larXlv:1101.5947l [hep-ph] 
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(1+ 2+, 3+, 4±) = -i|!/lJM^74 (m2, 0, 0, m^ s,,, S2s,m\ 0, 0, 0) 
{2sMmm) - + 2)(2 + 3)|7?4)) [32]m3 ^ 

(1+ 2+,3+,4t) {him^O,m')-l) 

m((7?i|(l + 2)(2 + 3)|r?4) + (r?ir?4)(2|l|2])[32] m((r?ir?4)(2|l|2] + (2??i)(3r?4)[32]) 

2(??il^)(??44^)(23)(2|l|2] 3(771 l^)(r?44^) (23)2 ' ^ ^ 



a[-^l(l+,2+,3+,4j^) = -i^4 (TO^0,0,m^Sl2,S23,0,TO^TO^m2) 

/ (2??i)(2?74) (2m^ + (2|1|2]) [32]2m3 (3??i)(3??4) (2m^ + (2|1|2]) [32]^m3 
V 2(r,il^)(r744^)(23)(2|l|3] + 2(r;il^)(,744^)(23)(3|l|2] 

_ (2 (2m2-S23) (?yir/4)+2(r/i|(l + 2)(2 + 3)|7?4)) [32]m3 \ 
2(771 l^)(r?44^) (23) ) 

rt 2 2 2 2 n^ (2'^^~*23) (m?74)[32]m3 n 2 n 2 2^ 

-J3(.23,m ,m ,m ,m (^a^)(r744^)(23)(2|l|2] -^3(.i2,0,m ,0,m ,m ) 

/ (2(7?ir74)(23)+ 4(2r;4)(3r?i))[32]m'^ (3771) (3774) (2|l|3][32]m3 (2r/i)(2774)(3|l|2][32]m3 
I (7711^) (7,44^) (23)2 + (77il^)(r,44^)(23)2(2|l|2] (77il^)(7744^)(23)2(2|l|2] 



, (277i)(27;4)(2|l|2][32]7r73 (3771) (3774) (2| 1 12] [32]7n3 \ . n 2 2 2x 

+ (77il^)(7744^)(23)2(2|l|3] " (77il^)(,44^)(23)2(3|l|2] j " ^3 (.23, 0, 0, m ,m ,m ) 

/ (2r/i)(2774)[32]2TO3 (3771) (3774) [32]2m3 (7717/4) [32]7773 \ 

1^2(77il^)(7/44^)(23)(2|l|3] " 2(771!^) (7/44^) (23) (3|1|2] (771!^) (7,44^) (23) j 

-F2(si2,0,77i2) 



2. (2si2(7/l7/4) - (77l|(l + 2)(2 + 3)|774)) [32]77l3 



(77il^)(r?44^)(23)(2|l|2]2 

+ UfHtXXAt) (12(^2,0,7772) - 1) + r.((77l|(l+ 2)(2 + 3)1774) + (771^^^^^^ 

2 ^ * ' ' ' * ^ 2V , : ) j 2(77il^)(7?44^)(23)(2|l|2] ' ^ ' 



-7A[^l(l+,2+,3+,4+)=- 



2m((77i774)(2|l|2] + (2771) (3774) [32]) 
(771 1^) (7744^) (23)3 [32] 



X I S23m%l3 (s23, 0, 0, m%, m%, m%) + 2m'jjF2 (s23, rnjj, m%) + ^.§23 j • (20) 
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